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Abstract: We examine whether future, nearly ah-sky galaxy redshift surveys, in combination 
with CMB priors, will be able to detect the signature of the cosmic neutrino background and 
determine the absolute neutrino mass scale. We also consider what constraints can be imposed 
on the effective number of neutrino species. In particular we consider two spectroscopic strategies 
in the near-IR, the so-called "slitless" and "multi-slit" approaches, whose examples are given by 
future space-based galaxy surveys, as EUCLID for the slitless case, or SPACE, JEDI, and possibly 
WFIRST in the future, for the multi-slit case. We find that, in combination with Planck, these 
galaxy probes will be able to detect at better than 3-sigma level and measure the mass of cosmic 
neutrinos: a) in a cosmology-independent way, if the sum of neutrino masses is above 0.1 eV; 
b) assuming spatial flatness and that dark energy is a cosmological constant, otherwise. We find 
that the sensitivity of such surveys is well suited to span the entire range of neutrino masses 
allowed by neutrino oscillation experiments, and to yield a clear detection of non-zero neutrino 
mass. The detection of the cosmic relic neutrino background with cosmological experiments will be 
a spectacular confirmation of our model for the early Universe and a window into one of the oldest 
relic components of our Universe. 
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1. Introduction 

Atmospheric and solar neutrino experiments have demonstrated that neutrinos have mass, 
implying a lower limit on the total neutrino mass given by M^, = Y^rriu ~ 0.05 eV |[l|. 
This is a clear indication that the standard model for particle physics is incomplete and 
that there must be new physics beyond it. The neutrino mass splitting required to explain 
observations of neutrino oscillations indicates that two hierarchies in the mass spectrum 
are possible: two light states and a heavy one (normal hierarchy, NH, with > 0.05 eV), 
or two heavy and one light (inverted hierarchy IH, with Mjy > 0.1 eV). A third possibility 
is that the absolute mass scale is much larger than the mass splittings and therefore the 
mass hierarchy does not matter (degenerate neutrino mass spectrum). 
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On-going and forthcoming neutrino experiments aim at determining the parameters of 
the neutrino mixing matrix and the nature of the neutrino mass (Dirac or Majorana). These 
experiments are sensitive to neutrino flavor and mixing angle, and to the absolute mass 
scale for large neutrino masses. As an example, beta-decay end-point spectra are sensitive 
to the neutrino mass, regardless of whether neutrinos are Dirac or Majaorana particles, 
and, the current limit on the effective electron neutrino mass is < 2.2 eV, coming from 
the Mainz and the Troitsk experiments, while KATRIN is expected to reach a sensitivity 
of ~ 0.2 eV ||2|, ^, ^. Near future neutrino oscillation data may resolve the neutrino mass 
hierarchy if one of the still unknown parameters, which relates flavor with mass states, is 
not too small. However, if the mixing angle is too small, oscillation data may be unable to 
solve this issue. 

On the other hand cosmological probes are blind to flavor but sensitive to the abso- 
lute mass scale even for small neutrino masses (see Fig.l). In fact, a thermal neutrino 
relic component in the Universe impacts both the expansion history and the growth of 
structure. Neutrinos with mass < 1 eV become non-relativistic after the epoch of re- 
combination probed by the CMB, and this mechanism allows massive neutrinos to alter 
the matter-radiation equality for a fixed Q^nh'^- Neutrino's radiation-like behaviour at 
early times changes the expansion rate, shifting the peak positions in the CMB angu- 
lar power spectrum, but this is somewhat degenerate with other cosmological parame- 
ters. WMAP7 alone constrains M,y < 1.3 eV ||5| and, thanks to improved sensitivity to 
polarisation and to the angular power spectrum damping tail, forecasts for the Planck 
satellite alone give Mi, ~ 0.2 — 0.4 eV, depending on the assumed cosmological model 
and fiducial neutrino mass (e.g., ^ and references therein). Massive neutrinos mod- 
ify structure formation on scales k > knr = 0.018(miy/leV)^/^Om^/i/Mpc, where /cnr 
is the wave- number corresponding to the Hubble horizon size at the epoch Znr, when a 
given neutrino species becomes non-relativistic. In particular, neutrinos free-stream and 
damp the galaxy power spectrum on scales k larger than the so called free-streaming scale 
k{s{z) = 0.82H{z)/{l + z)"^ {m,y / leV)hM-pc~^ ||l|, thereby modifying the shape of the mat- 
ter power spectrum in a redshift-dependent manner (see Fig. ^ and e.g. [10, 11, 12, ^]). 
Therefore, much more stringent constraints can be obtained by combining CMB data with 



large-scale structure (LSS) observations. Ref. [0, |l^ showed that present data-sets yield 



a robust upper limit of Mi, < 0.3 eV, almost ruling out the degenerate mass spectrum; this 
result was later confirmed by [16, 17 1. 

The forecasted sensitivity of future large-scale structure experiments, when combined 
with Planck CMB priors, indicate that cosmology should soon be able to detect signatures 
of the cosmic neutrino background and determine the sum of neutrino masses (e.g. [18, 19, 
m, 21 [ and references therein). Since cosmology is only weakly sensitive to the hierarchy 
|2^, a total neutrino mass determination from cosmology will be able to determine the 
hierarchy only if the underlying model is normal hierarchy and < 0.1 eV (see e.g. 
Fig. HI). A detection of the cosmic relic neutrino background (RNG) with cosmological 
experiments^ would be a spectacular confirmation of our model for the early Universe 
and a window into one of the oldest relic components of our Universe besides the one 
^Recall that neutrino experiments are not sensitive to relic neutrinos, as current generation of exper- 
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Figure 1: Constraints from neutrino oscillations (shaded regions) and from cosmology. In this 
parametrisation the sign of the mass splitting specifies the hierarchy. The red triangles show the 
fiducial models explored in this work and the light blue vertical bands our forecasted errors (see 
§5). For fiducial values below 0.1 eV a LCDM model must be assumed to obtain a detection 
with > 2-a statistical significance. For higher fiducial M^, we can marginalise over dark energy 
parameters and still obtain tight errors on M^. 



represented by the stochastic gravitational wave background. This consideration prompts 
us to examine whether future galaxy redshift surveys probing LSS will be able to detect the 
signature of the neutrino background and to determine the neutrino absolute mass scale. 

Beyond neutrino mass, cosmology is also sensitive to the number of neutrino species. 
In the standard model for particle physics there are three neutrinos; they decouple early 
in the cosmic history and then contribute to the relativistic energy density (i.e. as if they 
were radiation) with an effective number of neutrino species N^q = 3.046 (e.g. [|l|) until 
they become non-relativistic. Cosmology is sensitive to the physical energy density of 
relativistic particles, which include photons and neutrinos: Or = ^-y + N^s^u, where 
and flu are the energy density in photons and in one active neutrino species, respectively. 
CMB observations have constrained exquisitely well O^, thus constraints in fir can be used 
to study neutrino properties. Deviations from N(,s = 3.046 would indicate non-standard 

iments do not have sufficient energy resolution to cleanly pin down the signature of the RNG. Anyway, 
the beta-decay end-point spectrum is in principle also sensitive to the RNG, and this can be foreseen as a 
plausible perspective for future experiments only if neutrinos have masses of order eV, thus in the so called 
degenerate scheme for neutrino masses, which is still allowed by all present data, though slightly disfavored 
by cosmological observations p2|. 
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neutrino properties or additional effective relativistic species. Wliile the motivation for 
considering deviations from the standard model in the form of extra neutrino species has 
now disappeared |25, ^ 27], departures from the standard iVefj value could arise from decay 



of dark-matter particles 30, 31|, early quintessence [^], or more exotic models p3| . 

Relativistic particles affect the CMB and the matter power spectrum in two ways: 



a) through their anisotropic stress [34, g], and b) through their relativistic energy den- 
sity which alters the epoch of matter radiation equality. The ratio of CMB peak heights 
constrains matter-radiation equality yielding a degeneracy between iVefj and ilm^^- This 



degeneracy can be lifted by adding either cosmic expansion history data |35, ^] or 
adding the large-scale shape of the matter power spectrum: the power spectrum turnover 
scale is also related to matter-radiation equality given by the parameter T ~ (note 
the different scaling with h compared to the CMB constraint). LSS surveys can yield a 
measurement, at the same time, of both the cosmic expansion history (via the Baryon 
Acoustic Oscillations (BAO) signal), and the large scale turnover of the power spectrum. 
Present constraints are already competitive with nucleosynthesis constraints, and future 
data will offer the possibility to test consistency of the standard paradigm for the early 
Universe. In fact, nucleosynthesis constraints rely on physics describing the Universe when 
its energy scale was T ~ MeV, while cosmological constraints rely on physics at T ~ eV . 

In this paper we forecast errors on the total neutrino mass My and the effective number 
of relativistic species A'gff by combining Planck priors with data from future space-based 
galaxy redshift surveys in the near-IR. In particular, we consider two main survey strategies: 

• The first approach is to use "multi-slit" spectroscopy aimed at observing a pure 
magnitude-limited sample of galaxies selected in the near-IR (e.g. in the H-band at 
1.6 ^m) with a limiting magnitude appropriate to cover the desired redshift range. 
Examples of this approach are given by instruments where the efficient multi-slit 



capability is provided by micro-shutter arrays (MSA) (e.g. JEDI^ 39, |4^), or 
by digital micromirror devices (DMD) (e.g. SPACE |^ and possibly WFIRST^ in 
the future). With the multi-slit approach, all galaxy types (from passive ellipticals 
to starbursts) are observed, typically at < z < 2 — 3, if the observations are 
done in the near-IR, and provided that the targets are randomly selected from the 
magnitude-limited galaxy sample. 

The second approach is based on slitless spectroscopy (e.g. Euclid^ and JDEM^ 
P^, 43, 44 1) which, due to stronger sky background, is sensitive mostly to galaxies 



with emission lines (i.e. star- forming and AGN systems), and uses mainly Ha as a 
redshift tracer if the observations are done in the near-IR to cover the redshift range 
0.5 < z < 2. 



^http:/ /jedi. nhn.ou.edu/ 
^http:/ /wfirst. gsfc.nasa.gov/ 
''http:/ /sci.esa.int/euclid 
^http:/ /jdem. gsfc.nasa.gov/ 
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Forthcoming surveys will also have a weak gravitational lensing component, which will 
also be used to constrain neutrino properties (see e.g. |^). Here we concentrate on galaxy 
clustering as an independent and complementary probe. 

The rest of the paper is organised as follows. In § Q we review our method and the 
employed modelling. In § |3| we report the characteristics of the galaxy surveys considered 
in this work, and in § |^ we describe the adopted fiducial models and the explored space 
of cosmological parameters. In § |5| we present our results on the forecasted errors on the 
neutrino mass and number of neutrino species, and final in § ^ we draw our conclusions. 



2. Fisher matrix approach: P(A;)— method 



In this paper we adopt the Fisher matrix formalism to make predictions on neutrino masses 
and relativistic degrees of freedom from future galaxy redshift surveys. 

The Fisher matrix is defined as the second derivative of the natural logarithm of the 
likelihood surface about the maximum. In the approximation that the posterior distribution 
for the parameters is a multivariate Gaussian^ with mean fi = (x) and covariance matrix 



C = (xx*) 



/x/x , its elements are given by 



-Tr 



45|, |6 



47 



(2.1) 



where x is a N-dimensional vector representing the data set, whose components Xi are 
the fluctuations in the galaxy density relative to the mean in disjoint cells that cover 
the three-dimensional survey volume in a fine grid. The {9i\ denote the cosmological 
parameters within the assumed fiducial cosmology. 

In order to explore the cosmological parameter constraints from a given redshift survey, 
we need to specify the measurement uncertainties of the galaxy power spectrum. In general, 
the statistical error on the measurement of the galaxy power spectrum Pg{k) at a given 
wave-number bin is [Egjl 



P. 



2(27r) 



1 + 



1 



(2.2) 



' survey 

where Ug is the mean number density of galaxies, V^urvey is the comoving survey volume of 
the galaxy survey, and fi is the cosine of the angle between k and the line-of-sight direction 
fi = k ■ r/k. 

In general, the observed galaxy power spectrum is different from the true spectrum, and 
it can be reconstructed approximately assuming a reference cosmology (which we consider 
to be our fiducial cosmology) as (e.g. [pO[) 



-fobs(^i 



ref± ; ^ref 11 1 



Da{z)1,H{z) 



-fg(^ref±) ^ref II ) ^) + 



shot 



(2.3) 



^In practice, it can happen that the choice of parametrisation makes the posterior distribution slightly 
non-Gaussian. However, for the parametrisation chosen here, the error introduced by assuming Gaussianity 
in the posterior distribution can be considered as reasonably small, and therefore the Fisher matrix approach 
still holds as an excellent approximation for parameter forecasts. 



- 5 - 



where 



P. 



^ref ± ) ^ref II ; 



b{z) 



l + (3{z,k) 



Ay r I 

rci I 

'*^ref± + '^refll 



,(/c,z) 



(2.4) 



In Eq. (p.3|), -ff(-z) and -Da(-z) are the Hubble parameter and the angular diameter distance, 
respectively, and the prefactor {Da{z)'^^^H{z))/{Da{z)'^H{z)j.(,{) encapsulates the geomet- 
rical distortions due to the Alcock-Paczynski effect |5^. Their values in the reference 
cosmology are distinguished by the subscript 'ref, while those in the true cosmology have 
no subscript. k± and k\\ are the wave-numbers across and along the line of sight in the 
true cosmology, and they are related to the wave-numbers calculated assuming the refer- 
ence cosmology by k^ei± = k^DA{z) / DA{z)rei and k^^^ = k\\H{z)^cf/ H{z). Pshot is the 
unknown white shot noise that remains even after the conventional shot noise of inverse 
number density has been subtracted |^0[ , and which could arise from galaxy clustering bias 

), b{z) is the linear bias factor between 
fg{z, k)/b{z) is the linear redshift- 



even on large scales due to local bias |52]. In Eq. (| 
galaxy and matter density distributions, and (3{z, k) 
space distortion parameter [^], which in the presence of massive neutrinos depends on 
both redshift and wave-numbers, since in this case the linear growth rate fg{z,k) is scale 
dependent even at the linear level. We estimate fg{z, k) using the fitting formula of Ref. |p!T| 
(see the bottom-right panel of Fig. |2|). For the linear matter power spectrum -Pmatter(^5 z), 
we can encapsulate the effect of massive neutrino free-streaming into a redshift dependent 
total matter linear transfer function T{k,z) |54, 55, 56|, so that Pma.tter{k, z) in Eq. (^ 
takes the form 



Prr 



rik,z) 



87r2c^/coA^(A;o;^2 



T\k,z) 



Giz) 



G{z = 0) 



(2.5) 



where G{z) is the usual scale independent linear growth-factor in the absence of massive 
neutrino free-streaming, i.e. for k ^ (see Eq. (25) in Ref. [^6|), whose fiducial value in 
each redshift bin is computed through numerical integration of the differential equations 
governing the growth of linear perturbations in the presence of dark-energy ^7^. The 
redshift- dependent linear transfer function T{k, z) depends on matter, baryon and massive 
neutrino densities (neglecting dark-energy at early times), and is computed in each redshift 
bin using CAMB"^ |5^. As an example of its redshift dependence, in the top-left panel of 
Fig. ^ we consider the linear transfer function and show the ratio T{k, z)/T{k, z = 0) 
computed with CAMB at redshifts z = 0.5, 1, 1.5, 2 for a total neutrino mass M^, = 0.3 eV. 
On the other hand, in the top-right panel of Fig. ^, as an example of the neutrino free- 
streaming effect, we fix the redshift at z = and compute, for different neutrino masses, 
the ratio of the linear transfer function to the linear transfer function in absence of massive 
neutrinos. The power suppression due to neutrino free-streaming is evident and increases 
with the neutrino mass as well as the free-streaming scale. This suppression is also slightly 
dependent on the assumed mass hierarchy, as the blue-dotted and red-dashed lines clearly 
show. 



''http:/ /camb.info/ 
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Figure 2: Top Left: Ratio T{k, z)/T{k, z — 0) of the linear transfer functions computed with 
CAMB at redshifts z — 0, 0.5, 1, 1.5, 2, for a fiducial cosmology with a total neutrino mass = 0.3 
eV and a degenerate mass spectrum. Top Right: Ratio Tz=Q{k, Mi,) /Tz=o{k, Mi, = 0) between 
the linear transfer functions computed with CAMB for the different Af^-cosmologies described 
in Sec. 4 and the linear transfer function obtained assuming massless neutrinos. Bottom Left: 
Ratio Tz=o{k,Nes = 'i)/Tz=o{k,Nes = 3.03) between the linear transfer functions computed with 
CAMB assuming massless neutrinos and an effective number of relativistic species Ncs = 4 and 
A'cff — 3.04, respectively. Bottom Right: the function ii{k, f,y,ilde) = fg{M„ ^ 0)/fg{Mi, — 0), 
where — flit/flm- Note that fi{k, fi^,D,de) represents the scale dependent correction to fg{z), 
evaluated at M„ — 0.05. 



In Eq. ( ^.51) we have added the damping factor e ^'^i^'^'^l ^ due to redshift uncertainties, 



where 0"^ = {dr /dz)az, r{z) being the comoving distance [£9|, 50|, and we have assumed 



the power spectrum of primordial curvature perturbations, PTi{k), to be 

Al{k)^^-^^ = AUko)(^T , (2.6) 



where /cq = 0.002/Mpc, A|j(fco)|fid = 2.45 x 10^^ is the dimensionless amplitude of the 
primordial curvature perturbations evaluated at a pivot scale ko, and Ug is the scalar 



spectral index [gOl 



With the aim to make forecasts on the ability of future redshift galaxy surveys to 
constrain neutrino features, adopting different spectroscopic approaches (as discussed in 
§0), in this work we consider separately the effect on Pobs of the total neutrino mass Mjy 
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and the number of relativistic degrees of freedom A'eff, exploiting information from both 



the galaxy power spectrum shape and BAO distance indicators. In §5.6 we will analyse also 
the impact on neutrino mass constraints due to the inclusion of both growth-information 
and a Gaussian damping due to random peculiar velocities. 

In each redshift shell, with size Az = 0.1 and centred at redshift Zj, we choose the 
following set of parameters to describe Pohs{kref±, ^ref||) z): 

{H{zi),DA{zi),G(zi),P{zi,k),P'^f^^i,uJrn,(^b,C,ris,h} , (2.7) 

where C = ^eS ox C = = ^yh'^ (depending on the assumed fiducial cosmology, see fQ), 
Um = r^m^^; = f^fe/i^, where h is given by Hq = lOO/i km s~^ Mpc~^ , Hq being the Hub- 
ble constant. Vim, = M,y/(/i^93.8)eV, and 0.1, are respectively the total matter, massive 
neutrino, and baryon present-day energy densities, in units of the critical energy density 
of the Universe. Finally, since G{z), b{z), and the power spectrum normalisation Pq are 
completely degenerate, we have introduced the quantity G{zi) = {Po)^'^b{zi)G{zi)/G{zQ) 

In the limit where the survey volume is much larger than the scale of any features 
in Pobs(^); it has been shown that it is possible to redefine to be not the density 
fiuctuation in the n^^ spatial volume element, but the average power measured with the 
FKP method p| in a thin shell of radius kn in Fourier space. Under these assumptions 



the redshift survey Fisher matrix can be approximated as |46, 32 1 



TOO ainPobs(A?)ainPobs(^) .r. dk , , 

where the derivatives are evaluated at the parameter values pi of the fiducial model, and 
Ves is the effective volume of the survey: 



Vcsik,fi) 



ngPg{k,n) 



Purvey) (^'9) 



ngPg{k,fx) + 1 

where we have assumed that the comoving number density rig is constant in position. 
Due to azimuthal symmetry around the line of sight, the three-dimensional galaxy redshift 
power spectrum Pobs(^) depends only on k and /x, i.e. is reduced to two dimensions by 



symmetry |50]. 



To minimise nonlinear effects, we restrict wave-numbers to the quasi-linear regime, 
so that femax is given by requiring that the variance of matter fluctuations in a sphere of 
radius R is <j^(P) = 0.25 for R = 7r/(2A;max). This gives /cmax O.l/i Mpc~^ at z = 
and femax — 0.2/1 Mpc^^ at z = 1, well within the quasi-linear regime. In addition, we 
impose a uniform upper limit of /cmax < 0.2h Mpc~^ (i.e. femax = 0.2/i Mpc^^ at z > 1), to 
ensure that we are only considering the conservative linear regime, essentially unaffected 
by nonlinear effects. In each bin we adopt /cmin = 10~'^/i/Mpc, and we have verifled that 
changing the survey maximum scale kmin with the shell volume has almost no effect on the 
results. 
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For the moment, we do not include information from the amphtude G{zi) and the 
redshift space distortions I3{zi,k), so we marginaUse over these parameters^ and also over 
^Ihof Then we project p = {H{zi),DA{zi)^uJm-,oJh-,C-,f^s-,h} into the final sets q of cos- 
mological parameters described in §Q p3| , |64| . In this way we adopt the so-called "full 
P(A;)-method, marginalised over growth-information" |65|, and, to change from one set of 



parameters to another, we use |63| 



where F^^^ is the survey Fisher matrix for the set of parameters q, and -Fjj^^ is the survey 
Fisher matrix for the set of equivalent parameters p. 

We derive neutrino constraints with and without cosmic microwave background (CMB) 
priors; to this end we use the specifications of the Planck^ satellite. As explained in 
Appendix in order to describe CMB temperature and polarisation power spectra, we 
choose the parameter set 9 = {w^, t^;?,, C, 100^5, ln(10^'^A^(A:o)), n^, r}, where 0$ is the 
angular size of the sound horizon at last scattering, and r is the optical depth due to 
reionisation. After marginalisation over the optical depth, we propagate the Planck CMB 
Fisher matrix F^^^ into the final sets of parameters q, by using the appropriate Jacobian 
for the involved parameter transformation. 

The 1-0" error on qa marginalised over the other parameters is cr{qa) = \/ {F~^)aa, 
where F~^ is the inverse of the Fisher matrix. We then consider constraints in a two- 
parameter subspace, marginalising over the remaining parameters, in order to study the 
covariance between A'^gfr or My and the other cosmological parameters, respectively. 

Furthermore, to quantify the level of degeneracy between the different parameters, we 
estimate the so-called correlation coefficients, given by 



(2.11) 



where pa denotes one of the model parameters. When the coefficient |r| = 1, the two 
parameters are totally degenerate, while r = means they are uncorrelated. 

We evaluate (j{qa) and r both from survey data F^^^ , and from the combined 

7LSS I pCMB 
Q/3 "T ^ajS 



CMB+LSS data F^p = F^f + F^mb 



3. Surveys: two spectroscopic strategies 

In this work we forecast neutrino constraints using two different spectroscopic approaches. 
In particular, as mentioned in §||, we consider the cases relative to two space mission 
concepts under study: 

*In this case, we make derivatives of Poha{k,z) with respect to f3{zi,k). These derivatives are scale 
dependent, independently on the scale-dependence of /3. Then we integrate over k, as written Eq. (2.8); in 
this way we are left with the /3 redshift-dependence alone. Finally, we marginalise over it. 

^www.rssd.esa.int/index.php?project=planck 
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a EUCLID-like survey of Ha emission line galaxies, based on slitless spectroscopy of 
the sky. We adopt the empirical redshift distribution of Ha emission line galaxies 
derived by from observed Ha luminosity functions, and the bias function derived 
by using a galaxy formation simulation. In particular, we choose a flux limit of 
4x 10~^®ergs~^cm~^, a survey area of 20,000 deg^, a redshift success rate e = 0.5, a 
redshift accuracy of (Tz/{^ + z) < 0.001, and a redshift range 0.5 < z < 2.1. The total 
number of galaxies with redshift errors o"z/(l + z) < 0.001 from a slitless survey is 



well approximated by [35| 



% = 276.74 1^ ^ (/)-°-'(^")"' , (3.1) 
106 20000 0.5 ^ ^ ^ 

where / = //[10~^6ergs~^cm~^]^''. For this type of space-based slitless redshift 
survey we add in our forecasts also information from the ongoing Sloan Digital Sky 
Survey III (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS)^^ of luminous 
red galaxies (LRG). For this galaxy survey we assume that the LRG redshifts are 
measured over 0.1 < z < 0.5 with standard deviation ^^/(l + z) = 0.001, for a 
galaxy population with a fixed number density of n = 3 x 10^^/i^Mpc^^, and a fixed 
linear bias of 6 = 1.7 [^], over a survey area of 10,000 deg'^. 

a H-band magnitude limited survey of randomly sampled galaxies enabled by multi- 
slit spectroscopy (e.g., SPACE gl|], JEDI g], and possibly WFIRST in the 
future). To predict galaxy densities for such surveys we use the empirical galaxy 
redshift distribution compiled by Zamorani et al. from existing data l43], and we 



use predictions of galaxy bias from galaxy formation simulations |67]. We consider 



multi-slit surveys with a limiting magnitude of 77^^=22, a redshift success rate of 
90%, a sampling rate of 35%, a survey area of 20,000 deg^, a redshift accuracy of 
c^/(l + z) < 0.001, and a redshift range 0.1 < z < 2.1. The total number of galaxies 
with redshift errors ^"^/(l + z) < 0.001 from a multi-slit survey is well approximated 

by m 

^ = [192.21 + 197.03 {Hab - 22)^-^] . (3.2) 

106 L ^ \ AB J i 20000 0.9 X 0.35 ^ ^ 

Note that BOSS data are not added to the multi-slit galaxy redshift survey, since the 
latter has redshift ranges that extend to z ~ 0.1 [41, 43]. The case is different for Ha 



flux selected galaxies observed from space, since a wavelength range between 1 and 2 ^m 



naturally imposes a redshift range 0.52 < z < 2.05 in which Ha will be visible [43|. 

Furthermore, the bias functions for Ha flux and H-band magnitude selected galaxies 
increase with redshift, with the former being less strongly biased than the latter [B7|. In 



fact, the H-band traces massive structures (similar to selecting galaxies in the K-band), 



^°We note that this case is similar also to JDEM and ADEPT, 
^^http:/ /www. sdss3.org/surveys/boss.php 

^^For BOSS the actual redshift range is 0.1 < z < 0.7, but we do not take into account the shell 
0.5 < z < 0.7 in order to avoid overlapping in z between the space- and ground-based redshift surveys 
discussed in this work. 
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which makes them strongly biased. Star forming galaxies (which are selected by Ha flux), 
on the other hand, appear to avoid the cores of clusters and populate the filaments of the 
dark-matter structure, making them less biased than H-band galaxies [|67| ]. 

In conclusion, multi-slit surveys allow accurate redshift measurement for a larger num- 
ber of galaxies (and these galaxies are more biased tracers of large-scale structure than 
star-forming galaxies), and over a greater redshift range (extending to z ~ 0.1) than slit- 
less surveys. This can improve the constraints on neutrino masses, and, moreover, the 
ability to split a galaxy catalogue into red and blue galaxies could provide an important 
diagnostic test of potential systematic errors when measuring neutrino masses [^] . How- 
ever, multi-slit surveys have substantially stronger requirements in instrumentation and 
mission implementation |41|. 



4. Fiducial cosmologies: My and A'cfr 

The Fisher matrix approach propagates errors of galaxy power spectrum measurements 



Eq. (2.2) into errors of the cosmological parameters which characterise the underlying 
fiducial cosmology. According to the latest observations (e.g. [^] and refs. therein), we 
assume the fiducial cosmological model adopted in the Euclid Assessment Study Report 
||4^ with the exception that we normalise to the amplitude of the primordial curvature 
perturbations A^(/co) instead of a^: $7^ = 0.25, h = 0.7, A^(/co) = 2.45 x 10~^, Q;, = 
0.0445, Ug = 1^^. We consider neither primordial gravitational waves nor a scale dependent 
component of the scalar spectral index, and assume the matter energy density VLm to include 
the neutrino contribution when neutrinos are non-relativistic 

nrn = ^c + ^b + ^y , (4.1) 

Moreover, we also assume dark-energy to be a cosmic fluid described by a redshift 
dependent equation of state 

w,e{z) = ^ (4.2) 

Pde{z) 

where pde and pde represent respectively the pressure and energy density of the dark-energy 
fluid. 

This in turn yields a redshift dependent dark-energy density 



Pde{z) = PdeiO) exp 



1 + ^' . 



(4.3) 



for which we take a fiducial value ^l^e = 0.75. Finally, to compute our forecasts on dark- 



energy parameters, we adopt the widely used linear dark-energy equation of state [70, 71] 



Wde{a) =wo + {l- a)wa, (4.4) 



^^CMB constraints give ris ~ 0.96, anyway we have verified that adopting the latter value of changes 
neutrino mass constraints by ~5%, and the dark-energy equation of state parameter errors by ~10%, at 
the maximum. 
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where a = 1/(1 + z) is the scale factor normalised to unity at present, and we assume as 
fiducial values wq = —0.95 and Wa = 0, which lie well within the current 95% C.L. limits. 
In what follows, the dark-energy figure of merit (FoM) will be computed in terms of the 
conventional FoM for {wo,Wa) as proposed by DETF to compare dark-energy surveys. 
We will constraint the following final set of eight parameters 

q = {rim, ^de, h, AT^{ko),Qb, Wo, Wa,ns} , (4.5) 

which constitutes what we call our "base" parameters'"^, to which we add the C, parameter 
as we explain below'^. 

In what follows, we will consider six different fiducial models matching our "base" 
fiducial ACDM cosmology, in which we adopt the same fiducial values for the eight "base" 
parameters. Then we specify the assumed cosmological models: 

• The first model assumes a A'^efr^cosmology, where neutrinos are effectively massless 
but the the number of relativistic species ( = N^s can deviate from the standard 
value Ncs = 3.04. In this case the fiducial value A^efflfid = 3.04 is chosen, fixing 
= const = 0^6 H]. iVcff is given by the energy density associated to the total 
radiation 

n,. = n^{l + 0.22717Veff) , (4.6) 
where = 2.469 x lO^^/i^^ is the present-day photon energy density parameter for 



Te^b = 2.725 K PI 



The remaining models assume a M,.— cosmology, where A^cff is fixed at the fiducial 
value and C = is allowed to vary accordingly to the assumed fiducial neutrino 



mass spectrum [23]. In this case, we choose the following five fiducial values for the 



total neutrino mass consistent with current data ||7|, |5|: 

{0.3, 0.2 eV for degenerate spectrum 
0.125 eV for inverted hierarchy (4.7) 
0.125,0.05 eV for normal hierarchy 

The motivation for considering different fiducial models is that the dependence of the 
power spectrum on Mi, is nonlinear and thus the size of the forecasted error bar on M^, 
depends on the fiducial value chosen. 

At the CMB level, if neutrinos are still relativistic at the decoupling epoch, z ~ 1090, 
i.e. if the mass of the heaviest neutrino specie is tjIu < 0.58 eV, massive neutrinos do not 
affect the CMB power spectra, except through the gravitational lensing effect |7^, ^, |^, 
and, as a consequence, the dark-energy equation of state Wde is not degenerate with the 



^*Let us notice the adopted full P(fc)-method marginalised over growth-information does not give any 
constraint on A|j(fco), since the normalisation of the galaxy power spectrum is marginalised over. Therefore, 
the A Ij(fco) -errors shown in this work are forecasts from Planck alone. 



Finally, in §5.6 inchiding growth-information we will constrain the matter spectrum normalisation as 



together with the parameters in Eq. (4.5) 

different choice for the fiducial M^, would not affect the forecasted errors on A'off. 
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Table 1: Parameter 1-cr errors for slitless spectroscopy 



slitless+BOSS 



fiducial- 


> eV AU=0.2 eV 


M^=0.125 eV' M,.=0.125 eV" M^=0.05 eV' 


iVcff=3.04^ 




0.0140 


0.0137 


0.0139 


0.0138 


0.0137 


0.0124 




0.0260 


0.0265 


0.0258 


0.0257 


0.0253 


0.0256 




0.0032 


0.0031 


0.0032 


0.0032 


0.0032 


0.0034 


h 


0.0116 


0.0112 


0.0113 


0.0114 


0.0113 


0.0137 




0.1459 


0.1461 


0.1795 


0.1435 


0.1428 
















0.5435 


ns 


0.0233 


0.0225 


0.0314 


0.0228 


0.0220 


0.0326 


Wo 


0.0815 


0.0837 


0.0812 


0.0808 


0.0801 


0.0807 


Wa 


0.3461 


0.3573 


0.3450 


0.3440 


0.3386 


0.3324 


FoM 


51.75 


48.44 


52.44 


52.68 


54.25 


55.23 


slitless+BOSS+Planck 




0.0034 


0.0035 


0.0033 


0.0035 


0.0035 


0.0031 


^de 


0.0064 


0.0064 


0.0062 


0.0064 


0.0065 


0.0064 




0.0006 


0.0056 


0.0006 


0.0005 


0.0005 


0.0005 


h 


0.0043 


0.0043 


0.0043 


0.0043 


0.0042 


0.0046 




0.0347 


0.0433 


0.0311 


0.0441 


0.0526 
















0.0865 


ns 


0.0022 


0.0022 


0.0022 


0.0022 


0.0022 


0.0041 


Wo 


0.0732 


0.0721 


0.0716 


0.0715 


0.0709 


0.0705 




0.1760 


0.1742 


0.1713 


0.1725 


0.1722 


0.1664 




0.0250 


0.0226 


0.0226 


0.0227 


0.0244 


0.0227 


FoM 


245.07 


242.80 


259.32 


247.96 


240.31 


294.15 


'''for degenerate spectrum: 


mi ~ m2 ~ 


ms; ''for normal hierarchy: 


ms / 0, mi ? 


s m2 ~ 




'^for inverted hierarchy; mi 


^ m2, ma ; 


5i 0; '^fiducial cosmology with massless neutrinos 





neutrino mass. However, the limit on the the sum of neutrino masses degrades significantly 
when the dark-energy equation of state is a function of redshift as we assume in the present 
work, since dark-energy and massive neutrinos both affect the growth rate of structures 
1 75]. However, as we will show in the next section, the combination of CMB and LSS probes 
reduces or even breaks these degeneracies. The same does not happen for the number of 
relativistic species A'eff . Moreover, the degeneracies between Mi, and the other cosmological 
parameters can increase as the number of free parameters of the model increases, which 
could potentially bias the results for large k values As evident from Eqs. (pr^)-(p^), 
the model adopted in this paper falls within the two bias parameter models discussed in 
Ref. [^], which seem to mimic accurately the broad features of galaxy bias and redshift- 
space distortions from SDSS, leading to consistent constraints in the presence of massive 
neutrinos. 
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Table 2: Parameter 1-a errors for multi-slit spectroscopy 



multi-slit 



fiducial-^ 


M^==0.3 eV" 


AU=0.2 eV" 


Af^==0.125 eV' M, 


,=0.125 eV" 


AU=0.05 eV'' 


7Vcff=3.04'' 




0.0090 


0.0091 


0.0092 


0.0090 


0.0091 


0.0086 




0.0179 


0.0178 


0.0177 


0.0177 


0.0176 


0.0185 




0.0020 


0.0020 


0.0020 


0.0020 


0.0020 


0.0023 


h 


0.0079 


0.0078 


0.0079 


0.0079 


0.0079 


0.0098 




0.1229 


0.1113 


0.1321 


0.1110 


0.1126 
















0.4059 


ns 


0.0190 


0.0160 


0.0221 


0.0158 


0.0151 


0.0242 


Wo 


0.0617 


0.0614 


0.0619 


0.0613 


0.0612 


0.0629 


Wa 


0.2399 


0.2399 


0.2411 


0.2400 


0.2389 


0.2427 


FoM 


94.16 


94.78 


94.49 


95.12 


96.02 


94.09 


multi-slit -|-Planck 




0.0026 


0.0027 


0.0026 


0.0027 


0.0028 


0.0025 


^de 


0.0050 


0.0051 


0.0050 


0.0051 


0.0053 


0.0054 




0.0004 


0.0004 


0.0004 


0.0004 


0.0004 


0.0004 


h 


0.0033 


0.0033 


0.0033 


0.0033 


0.0033 


0.0038 




0.0296 


0.0376 


0.0268 


0.0388 


0.0463 
















0.0817 


ns 


0.0021 


0.0021 


0.0021 


0.0021 


0.0021 


0.0039 


Wo 


0.0554 


0.0552 


0.0552 


0.0552 


0.0551 


0.0551 


Wa 


0.1307 


0.1311 


0.1294 


0.1309 


0.1320 


0.1274 




0.0247 


0.0242 


0.0224 


0.0224 


0.0237 


0.0226 


FoM 


391.39 


376.55 


399.17 


377.41 


359.33 


441.78 


'''for degenerate spectrum: 


mi ~ m2 ~ ms 


; ''for normal hierarchy: 


ma / 0, mi ? 


s m2 ~ 




'^for inverted hierarchy; mi 


^ m2, ma 0; 


'^fiducial cosmology with massless neutrinos 





5. Results 

In this Section we present the predicted 1-a marginalised errors and correlations for the 
cosmological parameters considered in this work, focusing on the total neutrino mass M^, 
and the number of relativistic species A'efr, and comparing the results between the two spec- 
troscopic strategies described in We show forecasts from LSS alone and in combination 
with Planck priors. 

In Tables we show the marginalised errors for the six fiducial cosmologies considered 
here, computed adopting slitless and multi-slit spectroscopy, respectively; in Tables ^§ we 
report the corresponding correlation coefficients. Both LSS alone and LSS-I-CMB results 
are reported. 

Let us discuss the findings for the and A'^efj-cosmologies separately. 
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Table 3: Neutrino correlation coefficients for slitless spectroscopy 



slitless+BOSS 



fiducial- 


> eV M, 


.=0.2 eV 




=0.125 eV'' 


=0.125 eVM^ 


=0.05 eV'' 


iVcff=3.04' 




0.601 


0.573 




0.554 


0.549 


0.544 


-0.316 




-0.208 


-0.188 




-0.188 


-0.182 


-0.134 


0.221 




0.323 


0.331 




0.232 


0.303 


0.249 


-0.439 


h 


0.252 


0.259 




0.165 


0.230 


0.129 


0.602 




1 


1 




1 


1 


1 
























1 


ns 


0.567 


0.340 




0.717 


0.338 


0.364 


0.763 


Wq 


0.325 


0.310 




0.311 


0.298 


0.310 


-0.305 


Wa 


-0.463 


-0.437 




-0.424 


-0.417 


-0.413 


0.307 








slitlessH 


hBOSS+Planck 










0.344 


0.404 




0.331 


0.416 


0.485 


0.140 


^de 


-0.312 


-0.371 




-0.290 


-0.361 


-0.445 


0.389 




-0.141 


-0.131 




-0.136 


-0.117 


-0.169 


-0.062 


h 


0.019 


0.007 




0.023 


0.006 


0.019 


0.414 




1 


1 




1 


1 


1 




ns 


-0.089 


-0.071 




-0.034 


0.048 


0.020 


1 

0.844 


Wo 


0.011 


0.025 




0.008 


0.029 


0.031 


0.001 


Wa 


-0.166 


-0.202 




-0.157 


-0.204 


-0.266 


0.074 




0.106 


0.229 




0.073 


0.105 


0.288 


0.123 



'''for degenerate spectrum: mi « m2 ~ ms; ''for normal hierarchy: mg 7^ 0, mi ~ m2 ~ 
'^for inverted hierarchy: mi « m2, ms 0; ''fiducial cosmology with massless neutrinos 



5.1 cosmology: Correlations 

When considering forecasts from LSS alone, we find that M,^ is correlated with all the 
cosmological parameters affecting the galaxy power spectrum shape and BAO positions at 
scales k < kmax (see columns 2-6 in the upper panels of Tables In particular, there 

is a quite strong positive correlation r ~ 0.55 between the total neutrino mass My and 
the matter density 0^.. In fact, since neutrino free-streaming suppresses the total matter 
transfer function on scales smaller than the free-streaming scale kfg, increasing the neutrino 
mass produces on T{k, z) the opposite effect than increasing the total matter content of 
the Universe. For the same reason, massive neutrinos mimic the effect of a red tilt on 
the galaxy power spectrum, resulting in a positive correlation between M^^ and the scalar 
spectral index Ug. Moreover, Mj. is also positively correlated to the baryon density if 
the total matter content of the Universe is held fixed (a larger fif, enhances the BAO and 
the information content of T{k, z)). The correlation of M,^ with the remaining parameters 
h, Qde, Wq, and Wa can be explained looking at the Alcock-Paczynski prefactor in Eq. (^), 
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Table 4: Neutrino correlation coefficients for multi-slit spectroscopy 



multi-slit 



fiducial—)- 


AU=0.3 eV" 


A/^=0.2 eV" 


M^=0.125 eV' 


M^=0.125 eV" 


AU=0.05 eV'' 


iVeff=3.04'' 




0.607 


0.578 


0.577 


0.561 


0.544 


-0.447 




-0.151 


-0.144 


-0.136 


-0.136 


-0.134 


0.324 




0.284 


0.292 


0.217 


0.271 


0.249 


-0.551 


h 


0.175 


0.179 


0.110 


0.158 


0.129 


0.604 




1 


1 


1 


1 


1 
















1 


ns 


0.640 


0.442 


0.758 


0.434 


0.364 


0.814 


Wo 


0.341 


0.326 


0.342 


0.315 


0.310 


-0.368 


Wa 


-0.468 


-0.444 


-0.444 


-0.425 


-0.413 


0.407 


multi-slit+Planck 




0.355 


0.422 


0.337 


0.437 


0.485 


0.130 


^de 


-0.319 


-0.384 


-0.292 


-0.375 


-0.445 


0.472 




-0.201 


-0.184 


-0.195 


-0.167 


-0.169 


-0.049 


h 


0.052 


0.034 


0.058 


0.033 


0.019 


0.480 




1 


1 


1 


1 


1 




ns 


0.015 


0.013 


0.046 


0.029 


0.020 


1 

0.842 


Wo 


-0.011 


0.011 


-0.015 


0.016 


0.031 


0.005 




-0.168 


-0.217 


-0.156 


-0.222 


-0.266 


0.080 


Aliko) 


0.060 


0.177 


0.029 


0.068 


0.288 


0.153 



'^for degenerate spectrum: mi ~ m2 ~ ms; ''for normal hierarchy: mg 7^ 0, mi ~ 1x12 ~ 
'^for inverted hierarchy: mi « m2, ms 0; ''fiducial cosmology with massless neutrinos 



and recalling that for the moment we are not including information from the amplitude 
of the galaxy power spectrum, since we marginalise over the bias, the growth factor, the 
redshift space distortions, and the power spectrum normalisation. In this case we have 
to consider the expression of H{z) in presence of a non-vanishing spatial curvature of the 
Universe Ok = 1 — i^m — ^de ~ 

H{z) = Ho{^r[{l + zf - (1 + zf] + Vt^[{l + zf - (1 + zf] + 

Ode[(l + ^)3(l+"'0+«'a)g3^«4a-l) _ (1 + ^f]^ + (1 + zff^, (5.1) 

and remember that Da{z) is related to the inverse of H{zY'^ Given the fiducial values 
of /i, Orfe, wq, and Wa, and varying each parameter at a time, we see from Eq. ( ^?l| ) 

^^Let us specify that, at the redshifts covered by the LSS surveys considered in this work, we can safely 
consider that neutrinos are already non-realtivistic, so that they contribute to the Q.m entering in H{z). In 
fact, since we consider a minimum neutrino mass m^^i = 0.05 (normal hierarchy), the minimum redshift at 
which the massive neutrino component becomes non-relativistic is z,ir,i ~ 93.5 >> Zmax — 2.1 (see Eq. (A7) 
of Ref 0). 
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Figure 3: 2-paraineter M^-q^ joint contours with = ^de, wq, Wa, for the fiducial model with 
Ml, = 0.3 eV and a degenerate neutrino mass spectrum, obtained after combining the slitless survey 
data with BOSS data and Planck priors. The blue dotted line, the red dashed line and the orange 
dot-dashed line represent the 68% C.L., 95.4% C.L. and 99.73% C.L., respectively. The black solid 
line shows the 1-parameter confidence level at 1-a. 



that increasing h or wo enhances the observed power spectrum Pobs in contrast to the 
suppression induced by the increase of the total neutrino mass, so that these parameters 
are positively correlated with M^. On the contrary, from Eq. and looking at the sign 
of the term [(1 -|- 2)3(i+"'o+'i^<i)g3«)a(a-i) — (1 -|- z)'^], we also deduce that increasing Clde or 
Wa produces on Fobs the same effect as a larger neutrino mass, so these parameters are 
negatively correlated or "anti-correlated" to AIu. 

Except for the Ug case, we find that the level of correlation is on the average stable 
against the value of the fiducial total neutrino mass and the spectroscopic strategy adopted. 
For what concerns the mass hierarchy, at a given Mj^lgd, the Ug-My correlation looks to be 
larger by ^ 54% — 72% for normal hierarchy in comparison to the inverted one, at least when 
information from LSS alone are used, and the effect is more evident for slitless spectroscopy. 
In contrast, the h-My and Q^j-M^, correlations slightly decrease by ~ 34%— 24% respectively, 
for normal hierarchy compared to the inverted hierarchy. This can be understood if we 
consider that, for the same M^^^, in the case of normal hierarchy, the transfer function 
is slightly less suppressed on the scales of interest with respect to the inverted one (see 
Fig. 0, and § 5.2 for further comments on the relation between forecasted M,y-errors and 
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Figure 4: 2-paranieter Neff-qa joint contours with = ^de, WQ,Wa, ris for the fiducial model with 
extra relativistic degrees of freedom A^off = 3.04, obtained after combining the slitless survey data 
with BOSS data and Planck priors. The blue dotted line, the red dashed line and the orange dot- 
dashed line represent the 68% C.L., 95.4% C.L. and 99.73% C.L., respectively. The black solid line 
shows the 1-parameter confidence level at l-tr. 

the neutrino mass hierarchy). 

When Planck priors are added to the survey constraints, all degeneracies are either 
resolved or reduced, except for the covariance Miy-flde- In particular, the correlation be- 
tween Ml, and Ug is completely resolved, being reduced by ~ one order of magnitude. In 
some cases, the correlation coefficient r can even change sign (see columns 2-6 in the lower 
panels of Tables This change in the behaviour of r arises either due to the presence of 
dominant parameter degeneracies affecting the CMB spectrum, or because of marginalisa- 
tion of a high-dimension parameter space down to two variables. To summarise, after the 
inclusion of Planck priors, the remaining dominant correlations among Mi, and the other 
cosmological parameters are M^^-Qde, Afj/-^m, and M^-Wa- 

5.2 Mjy— cosmology: Forecasted error-bars 

The 1-0" errors of the parameters are shown in columns 2-6 of Tables We see that, 
with respect to the slitless spectroscopy, the multi-slit spectroscopy is able to reduce the 
neutrino mass errors of about 20%-30%, depending on the fiducial neutrino mass, if LSS 
data alone are used. In addition, for the same Mi,\^^, the \—a error on total neutrino mass 



-18- 




Figure 5: 1-parameter confidence levels at 1-a for and qa with qa = flde,wo,Wa,i^m for the 
fiducial model with M^, = 0.05 eV and a neutrino mass spectrum with normal hierarchy, obtained 
after combining the survey data with Planck priors. The blue solid line and the red dashed one 
represent the slitles+BOSS- and multi-slit-surveys cases respectively, as described in Sec. 3. 



for normal hierarchy is ~ 17% — 20% larger than for the inverted one. It looks like that the 
matter power spectrum is less able to give information on the total neutrino mass when the 
normal hierarchy is assumed as fiducial neutrino mass spectrum. This is similar to what 
found in Ref. ||2^ for the constraints on the neutrino mass hierarchy itself, when a normal 
hierarchy is assumed as the fiducial one. On the other hand, when CMB information are 
included, the Mjy-errors decrease by ~35% in favour of the normal hierarchy, at a given 
-^i^lfid- This difference arises from the changes in the free-streaming effect due to the 
assumed mass hierarchy, and is in agreement with the results in Ref. [^], which confirms 
that the expected errors on the neutrino masses depend not only on the sum of neutrino 
masses, but also on the order of the mass splitting between the neutrino mass states. 

When Planck priors are added, we find that the 1-a errors on M^y are in the range 
0.03 — 0.05 eV, depending on the fiducial total neutrino mass, with an average difference of 
15% between the two spectroscopic strategies, favouring again the multi-slit spectroscopy. 
This means that fixing some of the model free parameters, e.g. assuming a ACDM Universe, 
future spectroscopic galaxy surveys, combined with CMB probes, will be able to measure 
the minimum total neutrino mass M^, = 0.05 eV required by oscillation experiments; we 
will further comment on this in §5.6. Finally, depending on M;^|fid) the total CMB-I-LSS 
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Figure 6: 1-parameter confidence levels for N^s and qa with qa = ^lde,WQ,Wa,ns for the fiducial 
model with A'off = 3.04, obtained after combining the survey data with Planck priors. The blue solid 
line and the red dashed one represent the slitles+BOSS- and multi-slit-surveys cases respectively, 
as described in Sec. 3. 

dark-energy FoM decreases only by ~ 15% — 25% with respect to the value obtained if 
neutrinos are supposed to be massless, meaning that the "P(/c)-method marginalised 
over growth-information" is quite robust in constraining the dark-energy equation of state. 

5.3 A'cff cosmology: Correlations 

Likewise to the My case, we compute the 1-a errors and the correlation coefficients among 
A'eff and the cosmological parameters considered in Eq. (|4.5| ), for LSS alone and in combi- 
nation with Planck errors. 

Interpreting the sign of the correlations is not so straightforward for the A'eff- 
cosmology, since the number of relativistic species gives two opposite contributions to 
Pobs! and the total sign of the correlation depends on the dominant one, for each single 
cosmological parameter. In fact, from the bottom-left panel of Fig. ^, it is clear that a 
larger N^s value suppresses the transfer function T{k) on scales k < fcmax- On the other 
hand, from Eq. (|4.6| ) and Eq. ( |5.lD , we see that a larger A^eff value also increases the Alcock- 
Paczynski prefactor in Pobs- For what concerns the dark-energy parameters O^g) wq, Wa, 
and the dark-matter density ilm, we find that the Alcock-Paczynski prefactor dominates. 
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Figure 7: Top: M^-errors, for Afi^jfid = 0.05 eV, as functfons of the minimum redshift Zmin of 
the surveys, where we have fixed the maximum redshift z^ax = 2.1 for both the spectroscopic 
strategies. The lowest minimum redshifts considered are Zmin = 0.5 and Zmm = 0.1 for the shtless 
and muhi-sht spectroscopies, respectively. Bottom: M^-errors, for Afi,|fid = 0.05 eV, as functions of 
the maximum redshift Zmax of the surveys, where we have fixed the minimum redshifts Zmin — 0.5 
and Zmin — 0.1 for the shtless and muhi-sht spectroscopies, respectively. 



so that A'efT is positively correlated to Qde and Wa, and anti-correlated to ilm and wq. In 
contrast, for the other parameters, the T{k) suppression produces the larger effect and 
Nef[ results to be anti-correlated to Clb, and positively correlated to h and n^. See the 
last column in the upper panels of Tables for the LSS alone case. The degree of the 
correlation r is stable against the spectroscopic strategy adopted and is very large in the 
Ug-Ngfi case, being r ~ 0.8 with and without Planck priors. For the remaining cosmological 
parameters, all the correlations are reduced when CMB information are added, except for 
the covariance Ngff-Qdf,, as happens also for the M^z-cosmology. To summarise, after the 
inclusion of Planck priors, the remaining dominant correlations among A'efr and the other 
cosmological parameters are N^s-ns, Nes-Qde, and N^s-h. 

5.4 A'cflf cosmology: Forecasted errors 

The 1-0" errors of the parameters are shown in the last column of Tables [l[]2|. Also in this 
case, compared to the slitless spectroscopy, the multi-slit spectroscopy is able to reduce 
the A'^eff errors by '"^30% when LSS alone are used. When Planck priors are added, we find 
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Figure 8: Top: TVcff-errors as functions of the minimum redshift Zmin of the surveys, where we have 
fixed the maximum redshift Zmax =2.1 for both the spectroscopic strategies. The lowest minimum 
redshifts considered are Zmm = 0.5 and z„iin — 0.1 for the slitless and multi-slit spectroscopies, 
respectively. Bottom: A'cff-errors as functions of the maximum redshift z^ax of the surveys, where 
we have fixed the minimum redshifts Zmm — 0.5 and Zmin =0.1 for the slitless and multi-slit 
spectroscopies, respectively. 

a 1-a error on A'^efT of ~ 0.08, with a difference of only 6% between the two spectroscopy 
strategies, again in favour of the multi-slit one. 

5.5 Constraints on neutrino properties in the context of dark energy surveys 

The main science goal of the galaxy surveys considered in this work is to constrain dark 
energy. Considering neutrinos properties might degrade the dark energy constraints (be- 
cause of the introduction of extra parameters in the model) , or it may be that such surveys 
are not optimised (e.g., in their redshift coverage) to measure neutrino properties and thus 
perform sub-optimally for these parameters. 

We find that, depending on the fiducial Mi, value, the total CMB-I-LSS dark-energy 
FoM decreases only by ~ 15% — 25% with respect to the FoM obtained if neutrinos are 
assumed to be massless, and that, for the A'^efT~cosmology, the total CMB-I-LSS dark-energy 
FoM decreases only by ~ 5% with respect to the FoM obtained by holding N^s fixed. This 
means that the "P(/c)-method marginalised over growth-information" is quite robust in 
constraining the dark-energy equation of state. 
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In Figs. we show the jointly 2-parameter projected 68% C.L., 95.4% C.L. 
and 99.73% C.L. contours in the C,-qa sub-space, where C, = Mt^,Ncs and qa = 

Wq, Wa, f^jTt) Us, in the case of the shtless survey in combination with BOSS and Planck 
priors, for the A'eff- and M^(= 0.3 eV) -cosmologies. The black solid lines show the 1- 
parameter confidence levels at 1-a. The orientation of the ellipses reflects the correlations 
among the parameters shown in the lower panels of Tables 

Moreover, for a visual comparison of the constraints obtained with the two spectro- 
scopic strategies, in Figs. we show the 1-parameter confidence levels at 1-a in the (-Qa 
sub-space for the combinations slitless+BOSS+Planck and multi-slit -fPlanck, respectively. 

Finally, we consider the fiducial cosmology with Mi/|fid = 0.05 eV and, in the top 
panels of Fig. 0, we show the M,y-errors as functions of the of the 

surveys, with and without Planck priors. We fix the maximum redshift Zmax = 2.1 for both 
the spectroscopic strategies, while the lowest minimum redshifts are given by = 0.5 
and Zmin = 0.1 for the slitless and multi-slit spectroscopies, respectively. In the bottom 
panels of Fig. we show the M,y-errors as functions of the maximum redshift z^ax of the 
surveys, where we have fixed the minimum redshifts = 0.5 and = 0.1 for the 
slitless and multi-slit spectroscopies, respectively. For the slitless case, we have verified 
that extending the minimum redshift to Zmin = 0.1 changes the neutrino constraints by 
~0.8% only. We make the same analysis for the A'eff-errors and the corresponding trends 
are shown in Fig. ^. 

Note that concerning the possible degeneracies with the dark-energy parameters, we 
find that the major covariance is with ^de^ rather than with the dark energy equation of 
state parameters wq and Wa , and that extending the redshift range of the surveys considered 
would not reduce drastically the forecasted errors^^. 

5.6 The effects of growth inclusion and random peculiar velocities on neutrino 
mass constraints 

It is well known that galaxy peculiar velocities produce redshift-space distortions (RDS), 
which can be exploited by a large deep redshift survey to measure the growth rate of 
density fluctuations /g, within the same redshift bins in which H[z) is estimated via BAO. 
In particular, RDS allow to constrain fg times the normalisation of the power spectrum, 
i.e. /gCTg. In order to include growth information in our Fisher matrix analysis, following 



Ref. j65], we rewrite Eq. ( |2.3D as 

Da{z)1,H{z) 



-fofes(^ref±; ^ref 11 ) 



DA{zfH{zU 



<^fig{z) + fg{z,k)as{z) 



''^refll 



^ref± + ^refll 



X C{k,z) + Pshot, 
(5.2) 



where C{k,z) = Praa.ucrik, z)/ag{z), and (7sg{z) = b{z)as{z). We refer to this as the "full 
P(A;)-method, with growth-information included". Let us stress that we model RDS as an 

^^Concerning the dependence of M^- and A'cff-errors on the survey area, we note that F^^^ is Unearly 
dependent on the effective survey volume VcS (see Eqs. (2.8)-(2.9)), therefore a(Mv) and (t{Ncs), extracted 
from LSS data alone, are inversely proportional to the square root of the survey area, for a fixed redshift 
range. 
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additive component, using as free parameters {bag, fgds). fgag can be measured without 
knowing the bias b or the amphtude of the matter fluctuations as', therefore this parameter 
choice reduces possible systematic errors due to estimates of the bias |76|. 

While, for large separations, galaxy peculiar velocities give information on the growth 
of structures, on small scales random peculiar velocities cause the so-called Fingers of God 
(FoG), stretching compact structures along the line-of-sight |76, 77 1. Although, on the 
scales of interest in this work, this effect is expected to be moderate ||7^, |79|, we now include 
it in the Fisher analysis by introducing a Gaussian distribution for the pairwise velocity 
dispersion in configuration space, which produces a Gaussian damping ^-^^t^'^^l of the 
observed galaxy power spectrum Pobs- This effect is degenerate with possible inaccuracies 
in the observed redshifts due to a line-of-sight smearing of the structures, so we absorb 
it in the Gaussian damping factor of Eq. (|2.5| ). 

For the sake of simplicity, we consider only the fiducial cosmology with My = 0.05 
eV. In fact, in this case, the scale dependence of the growth rate fg due to free-streaming 
massive neutrinos can be assumed to be negligible. This is evident from the bottom- 
right panel of Fig. |2|, where we show the scale dependence correction to fg for My = 0.05 
eV, given by the function fi{k,fy,Qde) = fg{My / 0)/fg{My = 0), where fy = fly/Qm, 
introduced by Ref. [ [Tl] ] in their Eqs. (16)-(17). In this case the growth suppression is only 
of the order of 0.2%, affecting scales k > 0.1 /i/Mpc. 

Under these assumptions, when we include the effect of FoG in the Fisher analysis, 
we consider also (Jv{zi) as a scale independent variable, which is treated as a nuisance 
parameter to be marginalised over in each redshift bin, together with {asgizi), -Pg/joj}, 
and we project the errors on {H{zi), Da{z-i), fg{zi)as{zi), cOm, '^h-, C) f^s, h} into the 
final set of cosmological parameters {Vim-, ^de-, ^b-, wq, Wa, rig, ag}. This parameter 
set differs from Eq. (4^) only for the substitution A^(/co) o"8, whose fiducial value 
is fixed by A^(A;o) = 2.45 x 10~^ and the other parameters of the Mi,|fid(= O.OSeV)- 
cosmology described in ^ . We refer to this as the "full P(/c)-method, with FoG and 
growth-information included" . 

In what follows we compare the My constraints obtained by the three P(/i;)-methods 
considered in this work, for the fiducial cosmology with Mi/lgd = 0.05 eV. 



5.6.1 Growth inclusion effects 

We find that neutrino mass errors are quite stable at a{My) = 0.05 eV, against the adopted 
method (whether growth-information are included or marginalised over), and decrease only 
by 10%-20% when fgas measurements are included. We can understand this result as fol- 
lows. Qy affects the shape of the power spectrum, i.e. enters the transfer function T{k,z), 
which is sampled on a very large range of scales, including the P{k) turnover scale, by the 
nearly full-sky surveys under consideration. On the other hand, fg is only slightly depen- 
dent on My for A-Iy\fid = 0.05 eV (see the bottom-right panel of Fig. §). Consequently, the 
effect on the observed power spectrum shape dominates over the information extracted by 
measurements of fgcrs- This quantity, in turn, generates new correlations with My via the 
(Ts-term, which now we constrain simultaneously with the other cosmological parameters. 
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and which actually is anti-correlated with M,^^^ [ pO|] . 

On the other hand, if we suppose that early dark-energy is negligible, the dark-energy 
parameters ^de^ 'Wo and Wa do not enter the transfer function, and consequently growth 
information have relatively more weight when added to constraints from H{z) and D^iz) 
alone. 

As a result, the Mj^-errors are quite insensitive to growth inclusion, hence almost 
independent of the adopted P(A;)-method. This is in contrast to the dark-energy parameter 
constraints (see e.g. Ref. |11|). In fact, we find that, with respect to the "full P(fc)-method, 



marginalised over growth-information", the "full P(fc)-method, with growth-information 
included" is able to increase the dark-energy FoM by ~ 50% and ~ 60% for the slitless 
and multi-slit strategies, respectively, from both LSS data alone and in combination with 
Planck priors, when massive neutrinos are assumed in the fiducial cosmology. 

To summarise, we find that, due to the slight dependence of fg on Mi, when M,y|fid = 
0.05 eV, and due to the further degeneracy with ag (correlated also with the dark-energy 
parameters), we do not find a total effective gain on the accuracy of Mi, measurements, 
when growth-information are added. On the other hand, the value of the dark-energy 
FoM does increase when growth-information are included, even if it decreases by a factor 
~ 2 — 3 with respect to cosmologies where neutrinos are assumed to be massless, due to the 
correlation among Mjy and the dark-energy parameters. As confirmation of this degeneracy, 
we find that, when growth-information are added and if the dark-energy parameters ^de, 
wq, Wa are held fixed to their fiducial values, the errors a{My) decrease to 0.027 eV and 
0.025 eV, for the slitless and multi-slit spectroscopies combined with Planck, respectively. 

5.6.2 Incoherent velocity inclusion effects 

We expect that dark-energy parameter errors are somewhat sensitive to FoG effects. This 
can be understood in terms of correlation functions in the redshift-space; the stretching 
effect due to random peculiar velocities contrasts the flattening effect due to large-scale 
bulk velocities. Consequently, these two competing effects act along opposite directions on 
the dark-energy parameter constraints. 

We find that the dark-energy FoM obtained with the "full P(A;)-method, with FoG and 
growth-information included" result to be FoM^^^ = 51, 95 and FoM^^^+^'^^ = 268, 391, 
for the slitless and multi-slit spectroscopies respectively, i.e. very similar to the FoMs 
obtained from the "full P(A;)-method, marginalised over growth-information", as shown 
in Tables |l[^ which therefore can be considered a more stable approach against galaxy 
peculiar velocity uncertainties. 

On the other hand, the neutrino mass errors are expected to be almost stable at 
a{Mi,) = 0.05 when FoGs effects are taken into account by marginalising over cry(z), 
increasing only by 10%-14% with respect to the "full P{k)-met}iod, marginalised over 



^^In particular, using the "full P(fc)-method, with growth-information included", we find 
a{as) - 0.011,0.0085 for BOSS+slitless and multi-slit, respectively, and aiag) ~ 0.0014,0.0013 for 
BOSS+slitless+Planck and multi-slit+Planck, respectively. Using the "full P(fe)-method, with FoG and 
growth-information included", we find a{as) ~ 0.012, 0.0096 for BOSS+slitless and multi-slit, respectively, 
and again a{as) ~ 0.0014, 0.0013 for BOSS-|-slitless+Planck and multi-slit+Planck, respectively. 
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growth-information". Moreover, in this case, if the dark-energy parameters i^de^ wq, Wa 
are held fixed to their fiducial values, the errors a{Mjj) become 0.029-0.028 eV, for the 
slitless and multi-slit spectroscopies combined with Planck, respectively. In other words, 
we get errors which are only '^11% larger than the ones obtained without marginalising 



over (Jy{z) under the same assumptions on ^de-, wq, and Wa, as described in § 5.6.1 . 
5.7 Systematic effects, non-linearities and bias 

Fisher-matrix based forecasts are not particularly well suited to quantify systematic effects. 
The errors reported so far are statistical errors, which are meaningful only as long as 
they dominate over systematic errors. It is therefore important to consider sources of 
systematics and their possible effects on the recovered parameters. Possible sources of 
systematic errors of major concern are the effect of non-linearities and the effects of galaxy 
bias. In our analysis so far we have used the linear theory matter power spectrum and 
applied scale-independent bias to it. 

The description of non-linearities in the matter power spectrum in the presence of 
massive neutrinos is a relatively new subject. It has been addressed in several different 
ways: Refs. 75, |8^, use perturbation theory, Ref. [Q] used the time-RG flow 



approach and Refs. |85, |8g, g8|, 13| used different schemes of N-body simulations. From 
the above references it is clear that the effect of massive neutrinos on the matter power 
spectrum in the non-linear regime must be explored via N-body simulations to encompass 
all the relevant effects. Different simulations schemes and approximations agree already 
at or below the % level (for neutrino masses allowed by current observations and k < 1 
Mpc//i) indicating that non-linear effects on the matter power spectrum can in the future 
be modelled to the required accuracy. 

On the scales considered in this work the effects of non-linearities are small and 
statistical-errors could be further reduced considering smaller scales. As shown in |f75| , 
the effect of non linearities does not erase or reduce the effect of massive neutrinos. In fact, 
by comparing Fig. 4 of |^5[ and e.g., Fig. 4 of |^6[, it is apparent that the difference between 
the massless and massive neutrino case is enhanced by non-linearities. Thus while it will 
be mandatory to include non-linearities in the actual data analysis, the forecasted errors 
are not made artificially smaller by using the linear matter power spectrum to compute 
our Fisher matrices. 

Pushing to smaller scales however would worsen the systematic effect of scale- 
dependent and/or non-linear bias. Here we have made the simplifying assumption that 
bias is scale-independent up to kmax but the redshift dependence is not known and is 
marginalised over. Bias is known to be scale-independent on large, linear scales but to 
become non-linear and scale-dependent for small scales and/or for very massive halos. A 
scale-dependence of bias may mimic in part the effect of massive neutrinos. The scale- 
dependence of bias however is expected not to have the same redshift dependence as 
massive neutrino effects, thus offering the possibility to break a possible degeneracy. A 
scale-dependence of bias may cancel in part the effect of massive neutrinos. It arises be- 
cause halos, especially massive or rare ones, are non-linearly biased with respect to the dark 
matter. Scale-dependent bias (at least for dark matter halos hosting galaxies) is found first 
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Table 5: a{M^) and iT(A'off ) marginalised errors from LSS+CMB 









General cosmoloj 


lY 






fiducial^- M„ 


=0.3 eVM^ 


=0.2 eV" 


M^=0.125 eV'' 


M^=0.125 eV'M^ 


=0.05 eV'TV 


3ff=3.04'' 


slitless+BOSS+Planck 
multi-slit+Planck 


0.035 
0.030 


0.043 
0.038 


0.031 
0.027 


0.044 
0.039 


0.053 
0.046 


0.086 
0.082 








ACDM cosmolog 


y 






slitless+BOSS+Planck 
multi-slit+Planck 


0.017 
0.015 


0.019 
0.016 


0.017 
0.014 


0.021 
0.018 


0.021 
0.018 


0.023 
0.019 



"for degenerate spectrum: mi ~ m2 ~ ma; ''for normal hierarchy: ma 7^ 0, mi ~ m2 ~ 
'^for inverted hierarchy: mi « m2, ma ~ 0; '^fiducial cosmology with massless neutrinos 



to increase with increasing \k\ then decrease, but the bias of the galaxies hosted in the 
dark halos may be more complicated. This will be an important limitation in any practical 
application especially if we want to include mildly non-linear scales. There are, however, 
several possibilities to control or quantify systematics introduced by bias. In fact, the bias 
behaviour varies for differently selected objects (different colour or different brightness): 
splitting the sample in differently-biased tracers will thus help disentangle the systematic 
effect from the cosmological signal (approach similar to that of e.g.,[|69|]). Finally on large, 
linear scales, the neutrino mass splitting leaves a specific signature on the shape of the 
power spectrum that can also be used as a cross check of the My signal as illustrated in 
Ref. |87|. A more quantitative investigation of scale-dependent bias is beyond the scope of 
this paper, but we plan to study this important issue as more data and simulations become 
available. 



6. Conclusions 

In this work we have forecasted errors on the total neutrino mass My and the effective 
number of relativistic species N^.^, by combining Planck priors with future data from space- 
based spectroscopic galaxy redshift surveys in the near-IR. We have considered two survey 
strategies based on slitless and multi-slit spectroscopies. The assumed set of cosmological 
parameters is very general and takes into account a time-varying dark-energy equation of 
state, as well as a non- vanishing spatial curvature of the Universe. We exploited information 
from the galaxy power spectrum shape and BAO positions, marginalising over galaxy bias; 
thus our findings do not depend on bias measurement accuracy (as long as, on the large 
scales considered, bias is scale independent or its scale dependence is known), or modelling 
of the redshift dependence of bias [^] . 

The 1-0" errors are shown in Tables [l[]^, and the correlation coefficients in Tables |^^. 
In Figs. we show the joint 2-parameter confidence levels. 

Regarding M^^-errors, we find that the multi-slit spectroscopy is able to reduce the 
neutrino mass errors of about 20%-30% compared to the slitless spectroscopy, depending 
on the fiducial total neutrino mass, if LSS data alone are used. When Planck priors are 



-27- 



added, the 1-a errors on Mj^ are in the range 0.03 — 0.05 eV, depending on the fiducial 
neutrino mass, with an average difference of 15% between the two spectroscopic strategies, 
favouring the multi-sht spectroscopy. 

Moreover, depending on the fiducial Mj^-value, the total CMB+LSS dark-energy FoM, 
with growth-information marginalised over, decreases only by ~ 15% — 25% with respect to 
the value obtained if neutrinos are assumed to be massless (or their mass is assumed to be 
perfectly known), meaning that the "P(/c)-method marginalised over growth-information" 
is quite robust to assumptions about model cosmology when constraining the dark-energy 
equation of state. The situation is different when we include growth-information, since in 
this case the value of the dark-energy FoM decreases by a factor ~ 2 — 3 with respect to 
cosmologies that assume massless neutrinos. 

Considering the fiducial cosmology with Mi^lf^^ = 0.05 eV, in §^ we checked the 
stability of M,y-errors to the inclusion of growth-information and peculiar velocity uncer- 
tainties. We compared the following approaches: the "full P(A;)-method, marginalised over 
growth-information" , the "full P(A;)-method, with growth-information included" , and "full 
P(A;)-method, with FoG and growth-information included". We found that Mi,-errors are 
quite stable at a{My) = 0.05 eV, against the adopted method. This result is as expected, 
if we consider that, unlike dark energy parameters, My affects the shape of the power spec- 
trum via a redshift-dependent transfer function T(k,z), which is sampled on a very large 
range of scales including the P{k) turnover scale, therefore this effect dominates over the 
information extracted from measurements of /gcrg. 

Regarding A^eff~errors, again we find that, compared to the slitless spectroscopy, the 
multi-slit spectroscopy is able to reduce the A'gff-errors by ~30% when LSS alone are used. 
When Planck priors are added, we find a{Nc,s) ~ 0.08, with only a 6% difference between 
the two spectroscopy strategies, again in favour of the multi-slit one. The total CMB+LSS 
dark-energy FoM decreases only by ~ 5% with respect to the value obtained holding A'^eff 
fixed at its fiducial value, meaning that also in this case the "P(A:)-method marginalised 
over growth-information" is not too sensitive to assumptions about model cosmology when 
constraining the dark-energy equation of state. 

Finally, in Table | we summarise the dependence of the M^- and A'efr-errors on the 
model cosmology, for the two spectroscopic strategies combined with Planck. We conclude 
that, if My is > 0.1 eV, these surveys will be able to determine the neutrino mass scale 
independently of the model cosmology assumed. If My is < 0.1 eV, the sum of neutrino 
masses, and in particular the minimum neutrino mass required by neutrino oscillations, 
can be measured in the context of a ACDM model. 

This means that future spectroscopic galaxy surveys, such as Euclid or SPACE, JEDI, 
and possibly WFIRST in the future, will be able to cover the entire parameter space for 
neutrino mass allowed by oscillations experiments 

Moreover, as summarised in Fig. 0, they will be competitive with future 3D cosmic 
shear photometric surveys, which, in combination with Planck priors, will give similar 
constraints on My and N^s [^. Since, these two kinds of LSS probe are affected by different 
systematics, their constraints on neutrino masses and relativistic degrees of freedom will 
provide a consistency check of the two independent measurement methods. 
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We conclude that future nearly all-sky spectroscopic galaxy surveys will detect the 
cosmic neutrino background at high statistical significance, and provide a measurement 
of the neutrino mass scale. This will provide an important confirmation of our model for 
the early Universe, and crucial insights into neutrino properties, highly complementary to 
future particle physics experiments. 
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A. Planck priors 

In this work we use the Planck mission parameter constraints as CMB priors, by estimating 
the cosmological parameter errors via measurements of the temperature and polarisation 
power spectra. As CMB anisotropics, with the exception of the integrated Sachs- Wolfe 
effect, are not able to constrain the equation of state of dark-energy (wo^Wa)'^^ , we follow 



the prescription laid out by DETF [72] 



We do not include any B-mode in our forecasts and assume no tensor mode contribu- 
tion to the power spectra. We use the 100 GHz, 143 GHz, and 217 GHz channels as science 
channels. These channels have a beam of Sfwhm = 9.5', Sfwhm = 7.1', and 0fwhm = 5', 
respectively, and sensitivities of fj-r = 2.5fiK/K, ax = 2.2fiK/ K , ctt = A.^^K/K for 
temperature, and ap = AfiK/K, ap = A.2^K/K, ap = 9.8fj,K/K for polarisation, respec- 
tively. We take /sky = 0.80 as the sky fraction in order to account for galactic obstruction, 
and use a minimum i-mode imin = 30 in order to avoid problems with polarisation fore- 
grounds and not to include information from the late Integrated Sachs- Wolfe effect, which 
depends on the specific dark-energy model. We discard temperature and polarisation data 
at ^ > 2000 to reduce sensitivity to contributions from patchy reionisation and point source 



contamination (see [72| and references therein). 

We assume a ACDM fiducial cosmology, and choose the following set of parame- 
ters to describe the temperature and polarisation power spectra 9 = (a;m,a;fe,C,100 x 
^s, ln(10^''A^(A;o)), ng, r), where 9s is the angular size of the sound horizon at last scat- 
tering, and r is the optical depth due to reionisation. Note that a different parameter set 
is assumed in [39|. 



The Fisher matrix for CMB power spectrum is given by [90, 



I X,Y * ^ 



^"On the contrary, using {wo,Wa) as model parameters to compute the CMB Fisher matrix could artifi- 
cially break exiting degeneracies. 
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where 6i are the parameters to constrain, Cx,i is the harmonic power spectrum for the 
temperature-temperature {X = TT), temperature-E-polarisation {X = TE) and the E- 
polarisation-E-polarisation {X = EE) power spectrum. The covariance 

COVxY of the 

errors for the various power spectra is given by the fourth moment of the distribution, 
which under Gaussian assumptions is entirely given in terms of the Cx i with 



COVt.t 
CO\e,e 
COYte,te 

COYt,e 

C0\ T,TE 

COYe,te 



fi {Ct,i + W^^B^ 
fi{CE,i + Wp'Bf^)^ 



fi 



G 



TEI 



+ 



[Ct,i + W^^Bf"") {Ce,i + Wp^Bf^) 



TE,l 

hCTE,l {Ct,1 + ^^r^) 
hCTE,l {Ce,i + Wp^B^^) 



where fn 



(2^+l)/sky 



T,P"{wYiml 



(A.2) 
(A.3) 
(A.4) 

(A.5) 
(A.6) 
(A.7) 

being the weight per 



sohd angle for temperature and polarisation respectively, with a l-o" sensitivity per pixel 
of (T^p and a beam of extent, for each frequency channel c. The beam window 

function is given in terms of the full width half maximum (fwhm) beam width by Bf = 
Y.MfW-^p/WT,p. where {B<lf = exp , = ^(8 In 2) and 



/sky is the sky fraction [p2 |. 

We then calculate the Planck CMB Fisher matrix with the help of the publicly available 
C AMB code . Finally, we transform the Planck Fisher matrix for the DETF parameter 
set to the final parameter sets q considered in this work (see §Q and §^.6|), using the 
transformation 

pCMB _ pCMB 

dq^ Bn. ■ 



{A.i 
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